Size effects in dislocation depinning models for plastic yield 
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Typically, the plastic yield stress of a sample is determined from a stress-strain curve by defining 
a yield strain and reading off the stress required to attain it. However, it is not a priori clear that 
yield strengths of microscale samples measured this way should display the correct finite size scaling. 
Here we study plastic yield as a depinning transition of a 1 + 1 dimensional interface, and consider 
how finite size effects depend on the choice of yield strain, as well as the presence of hardening and 
the strength of elastic coupling. Our results indicate that in sufficiently large systems, the choice of 
yield strain is unimportant, but in smaller systems one must take care to avoid spurious effects. 
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I. INTRODUCTION 

A complete understanding of microscale plasticity is 
important both for fundamental and technological rea- 
sons, and much experimental and theoretical effort has 
been expended on this problem in recent years. A key 
result that has emerged — and that is not intuitively 
obvious from our everyday experience of classical bulk 
plasticity — is that on the microscale plastic flow is sub- 
ject to size effects, with recent experiments showing that 
in submicron samples^he yield stress increases with de- 
creasing sample size However, the origins of these 
finite size effects are not clear, and many mechanisms 
have been proposed (for reviews, see Refs 043)- 

In typical experimental studies, the yield stress is de- 
termined from measured stress-strain curves by defining 
a yield strain 7y, typically a few percent, and reading off 
the stress required to attain such a strain. This procedure 
suffers from the fact that measured stress-strain curves 
are not smooth and the transition to plastic yield is gen- 
erally preceded by power law distributed avalanches 
This raises the question of how the yield strain should be 
selected to ensure that the measured yield stress is not 
systematically over- or under-estimated in a way that 
affects finite size scaling. Furthermore, it is not a pri- 
ori clear that one can compare experimental results in 
which different threshold strains have been used, making 
it difficult to synthesize results from different studies. 

In this paper, we approach this problem from a theo- 
retical perspective, by studying numerical simulations of 
1 -|- 1 dimensional depinning at zero temperature. This is 
the simplest model for plastic yielding in a crystal: the in- 
terface represents a dislocation, driven through a random 
potential generated by disorder and other dislocations in 
the sample. Depinning has previously been studied in a 
variety of contexts, using both numerical and analytical 
techniques > and it is well-known that there exists a 

critical driving force /c, above which the interface is de- 
pinned and has finite velocity always, and below which 
the interface will always come to a stop at long times. 
This transition is analogous to plastic yielding, in which 
stresses above a finite threshold are able to induce large 



strain increases that ideally are limited only by sample 
failure. 

Like other phase transitions, the depinning transition 
is subject to finite size effects, so that the depinning force 
in a system of linear size L has the asymptotic form 



fc{L) ^foo + aL 



(1) 



where foo is the depinning force in the limit L ^ oo and v 
is the finite size exponent, determined by the universality 
class of the model. A key question that this work aims 
to address is whether the observed finite size scaling is 
affected by how the depinning transition is measured. 

In particular, we can imitate the experimental proce- 
dure for measuring yield stress by mapping our simula- 
tions of interface depinning onto plots of stress vs plastic 
strain from stress controlled experiments. The motion of 
a dislocation with Burgers vector b over distance 2; in a 
system of size results in a strain increase of 7 = bx/L^, 
and the driving force is proportional to an external stress. 
We can therefore consider a plot of the mean position of 
the interface vs driving force to represent plastic strain 
vs stress, as illustrated in Fig.[lja). In analogy to the ap- 
proach taken in experimental analysis, we define a target 
displacement A = •yL'^/b to define the onset of depinning, 
and study the effects of varying this target. In what fol- 
lows, we will measure lengths in units of b. 

We note that this approach to measure the depinning 
threshold raises another issue. It has previously been 
pointed out (igI . [itI . [l9j that in a simulation of a system 
of finite size L x M with periodic boundary conditions 
in both directions, a "natural" choice of M is M ~ L^, 
where ( is the interface roughness exponent and gives 
the interface width at criticality. In this way, one ensures 
that for all L, the system contains the same number of 
critical configurations. On the other hand, our target 
scales as L^, so although it is natural for the plasticity 
scenario at hand, the system aspect ratios increase with 
L. We will see that this affects the distribution of mea- 
sured depinning forces, as anticipated (l6l. [l7|. 

The depinning models can also be related to strain- 
controlled experiments, by driving the interface via a 
spring of stiffness k, and adjusting the externally-applied 
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stress so that the total strain — the elastic strain of the 
spring, partially relieved by the plastic deformation aris- 
ing from interface motion — increases smoothly. The 
set-up and a typical stress-strain curve are illustrated in 
Fig. [T{b). Again in this context, we study the effect of 
the choice of yield strain on measured yield stress, and 
discuss how these effects can be related to strain harden- 
ing in a stress controlled experiment. 

The outline of the paper is as follows. In Section |lTl 
we use the quenched Edwards- Wilkinson model for in- 
terface depinning as a case study for the role of tar- 
get choice on measured properties, considering not only 
the finite size scaling of the mean depinning force but 
also how the choice of target affects its whole distribu- 
tion. The ideas developed in Section are indepen- 
dent of universality class, which we demonstrate in Sec- 
tion [IlT] via results from the Kardar-Parisi-Zhang (KPZ) 
model and a mean field model. Finally, in Section IIVI 
we discuss two closely-related extensions to the quenched 
Edwards- Wilkinson model, namely strain hardening and 
the above-mentioned strain controlled experiments. We 
show that finite size effects in these models are more com- 
plex than in the other models we study, depending on 
both the finite size scaling of the depinning force and 
also the hardening amplitude or stiffness of the coupled 
spring. 

II. EFFECTS OF TARGET CHOICE: THE 
QUENCHED EDWARDS- WILKINSON 
EQUATION 

We start by considering a simple linear model for the 
interplay between elasticity, random pinning and an ex- 
ternal drive, namely, the quenched Edwards- Wilkinson 
model [1, 0| • This will allow us to discern the effect of 
the choice of target interface displacement in a straight- 
forward manner, giving results that will aid the interpre- 
tation of more complex models that follow. 

In general, an overdamped elastic interface character- 
ized by height function u{x) and subject to a random 
potential and a uniform external drive / is described by 
overdamped equations of motion 

iw(a:) = /tot (a;) (2a) 
= f + fci{x,u{x)) ~ T]{x,u{x)), (2b) 

where the damping constant F is hereafter set to unity. 
The self-interaction of the interface, /d, promotes it to be 
flat, and may be local or long-range, depending on the 
model at hand. The self-interaction is in competition 
with the quenched disorder ry, which causes the interface 
to prefer to pass through sites with low potential energy. 
In the quenched Edwards- Wilkinson model that we study 
in this Section, the elastic force at position x on the line 
is given by V^w(x). 

In this work, we study a cellular automaton version 
of this model, by discretizing x and u. The elastic force 
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FIG. 1: (a)Typical "stress-strain" (driving force, /, vs mean 
interface displacement, u) curves from simulations of the 
quenched Edwards- Wilkinson equation, with (circles) and 
without (triangles) hardening. The black solid line indicates 
the f (X u behavior of the model with hardening as m — ^ oo 
(see Section|lV|. Depinning forces are determined by selecting 
a target value of u and reading off the driving force required 
to attain this target, (b) Data from a hardening simulation 
recast as a strain controlled experiment with machine stiffness 
k, using the transformations ut = f — ku/L^, 7t = f /k, as 
described in Section HVl The yield stress, aT,y is determined 
by selecting a yield strain 7t,i/ and reading off the stress re- 
quired to attain it. 

is given by u^^i + u^-i — 2u^. We allow only forward 
motion of the interface: at each simulation step, a site 
x may increase Ux by one unit if the total force acting 
on it is positive, or remain stationary if the force is zero 
or negative. The random potential is implemented by 
drawing a random number for each x,Ux- In general, 
these are from a uniform distribution over [0, 1], but in 
Section IIIII we will also study Gaussian distributed dis- 
order, because when all sites are coupled to one another, 
the disorder distribution can be important We use 
periodic boundary conditions in the x direction, and open 
boundary conditions in the u direction. 

In order to generate "stress-strain" curves in a stress- 
controlled simulation, we start with a flat interface with 
Ux = everywhere and / = 0. We increase the drive / 
adiabatically, by flnding the site x* with the smallest re- 
quired driving force for motion, and increasing / so that 
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FIG. 2: (a) Mean measured depinning force (fc) plotted 
against system linear size L, from simulations of the quenched 
Edwards- Wilkinson equation. The target strains used are 
7 = 0.01 (circles), 0.07 (triangles), 0.25 (squares) and 2 (in- 
verted triangles). For all targets, fc approaches the same 
value as L ^ oo, but the behavior for smaller L depends on 
7. Error bars (typically smaller than the symbols) represent 
the standard error, (b) The same data, collapsed using Eq. ^ 
and (|3]). The black solid line is proportional to L""'^^, which 
is the expected finite size scaling for the model for large L. 

X* can move. This may trigger an avalanche via changes 
in /el, and we update sites in parallel until no more can 
move. We record the value of / and the mean position of 
the interface attained after it has stopped moving, and 
then increase the drive again. As mentioned in the Intro- 
duction, the record of driving forces and mean interface 
positions is the "stress-strain" curve, since the motion of 
a dislocation by distance x increases the plastic strain 
by x/L"^ (we work in units where the Burgers vector b is 
unity) . We continue to increase the drive until the mean 
interface position passes the target jL'^. The driving 
force required for the target to be attained is denoted fc, 
and for simplicity we will call it a "depinning" threshold, 
even if this is not strictly accurate. 

Figure [IJa) shows the mean force required for the in- 
terface to attain the target mean position 7^^, as a func- 
tion of system linear size L, for a range of 7 values be- 
tween 0.01 and 2. In the large-i limit, we find that 
fdL) approaches foo — 0.424 according to L"^/"^, with 
l/v = 0.75, in agreement with previous studies of finite 
size scaling [13: 13 • However, for small L, a distinct 7 
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FIG. 3: Standard deviation in the measured depinning force, 
a"(/c), plotted against system linear size L. For sufficiently 
large targets, a{fc) ~ L""'^^, as indicated by the solid black 
line. The target strains used are 7 = 0.07 (circles), 0.15 
(triangles), 2 (squares). Because the distributions of fc are 
not Gaussian (see Fig.[4l), the error bars have been estimated 
using the method in Ref. [2l| . 

dependence becomes apparent: for small 7, fc initially 
grows with L, but for large 7, fc{L) decreases monoton- 
ically. 

This result can be understood by considering that 
there are two length scales in the direction of interface 
propagation. The first is the interface width at depin- 
ning, which is given by wL'', with the roughness expo- 
nent C = 1-25 for the quenched Edwards- Wilkinson equa- 
tion ^ ^ constant. The second length scale is 
our choice of target, which is 7^^. The two length scales 
are equal when 



giving a crossover length scale Lc = {'w/^y^^^~''\ 

If L < Lc, the target is shorter than the interface 
width at criticality, and the simulation is halted during 
the early stages of interface motion, before the depinning 
transition can occur. In this regime, increasing L — and 
thereby the target — should increase the final driving 
force applied before the simulation is stopped. On the 
other hand, for L > Lc, we can expect the interface to 
undergo the depinning transition, and the standard finite 
size scaling of Eq. ^ should hold. At Lc, there is a tran- 
sition from the former behavior to the latter. However, 
if 7 is sufficiently large, the target is always larger than 
the interface width, that is, Lc — J' as 7 — 00. 

To verify this interpretation, we collapse the data from 
simulations with different 7 values. From Eq. ([3|), we see 
that we must rescale L — Lj~^/'-^~''\ Then from the 
finite size scaling ([T]) and the fact that = 1/(2 — Q , 
we rescale {fc — foo) by 7. As seen in Fig. [UJb), the 
data collapse is very good. Furthermore, that the correct 
finite size scaling laws are recovered for sufficiently large 
targets is further evidenced by the measured scaling of 
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FIG. 5: Data collapse for the KPZ model, based on the 
same arguments as the quenched Edwards- Wilkinson model 
(Eq.[l]and[3]), using fitted parameters C, = 0.635, foo = 0.226, 
1/v = 0.55. The symbols represent 7 = O.f (circles), 7 = 0.01 
(squares) and 7 = 0.005 (diamonds). The standard error is 
smaller than the symbol size. 



widths. In agreement with Refs. [la, ll7|, as 7 increases, 
the distribution approaches the Gumbel limit, albeit very 
slowly in the tails of the distributions. 



FIG. 4: (a) Histograms of the measured depinning force over 
10"" trials, for L = 100, with 7 = 0.01 (circles), 1 (squares), 10 
(diamonds) and 100 (triangles). The distributions have been 
normalized to have zero mean and unit variance. Error bars 
are smaller than the symbol size, (b) Cumulative distribution 
functions, showing convergence towards the Gumbel distribu- 
tion (a straight line on this double-log plot) as 7 increases. For 
comparison, in both plots, the Gaussian distribution (dashed 
line) and Gumbel distribution (solid thin line) are also shown. 



the standard deviation of fc{L), as plotted in Fig.[3l This 
has the asymptotic form L"^/"^, as expected from finite 
size scaling. 

We can go beyond the moments of the distribution of 
fc to consider the distribution itself. It has previously 
been shown [l^ [i3 that the distribution of fc measured 
in a finite sample of size L x M (with open boundary 
conditions in the M direction) depends on the ratio of 
M to the interface width wL''. Heuristically, the system 
can be divided into n = M/{wL'^) subsystems which are 
approximately independent, each with its own depinning 
force. The depinning force in the entire L x M system 
is the maximum of these, and if n — (X) the measured 
depinning force is described by extreme value statistics, 
in particular, the Gumbel distribution, which has cumu- 
lative distribution P{x) — exp(— e"'^"''-'/^-'), where fj, is 
the modal value of x and the standard deviation of the 
distribution is (3tt/^. In Fig. |4l we plot the distribu- 
tion of fc (rescaled to have zero mean and unit standard 
deviation) for L = 100, for a range of 7 values. The tar- 
gets range from a few interface widths to ~ IC* interface 



III. OTHER SIMPLE DEPINNING MODELS: 
KPZ AND MEAN FIELD 

The role of the parameter 7 relative to the interface 
width does not depend specifically on the choice of in- 
teraction /el used in the quenched Edwards- Wilkinson 
equation, and we have verified that models from other 
universality classes are affected by the choice of 7 in 
the same way. For example, in the KPZ model |23 |. 
where the elastic force fd = \7^u{x) is supplemented 
by a nonlinear term proportional to (Vm)^, we mea- 
sure a roughness exponent of C = 0.635 (consistent with 
other measurements [l^ [l^ ) and a finite size exponent 
of 1/1/ w 0.55. Per Eq. ©, fc{L) curves for different 
7 can be collapsed via the rescaling L — >■ 7^/*^^~^)L, 
{fc - foo) ^ 7"'/^"^^"'^^H/c - /oo), as shown in Fig.E] 

The quenched Edwards- Wilkinson and KPZ models 
both have local interactions, and the interface width 
scales with length as L"^ . We now turn to a mean- field 
model [13, [13, [13 , in which every site interacts with ev- 
ery other site, and the interface width does not depend 
on the length scale it is measured over [C, —Q). We cou- 
ple the interface to its center of mass, u. The total force 
on site x is 

fx = f + Jiu-Uj;) +riu,x- (4) 

In Fig. ini we show the finite size behavior of the de- 
pinning force for this model. The finite size exponent 
takes the mean- field value of v = 2. For fixed J we can 
collapse data from simulations with different 7 values fol- 
lowing the same argument as for the quenched Edwards- 
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FIG. 6: Rescaled depinning force vs target displacement 7L , 
in the mean field model. The solid black lines are proportional 
to L~^^^ . Closed symbols represent data for J = 1, open 
symbols represent data for J = 0.1, the 7 values are 0.005 
(squares), 0.01 (triangles), 0.1 (circles) and 1 (diamonds). For 
J = 0.1, /oo ~ 0.64, for J = 1, /cx, < 10"^ 



Wilkinson model, but with the interface width constant. 
As seen in Fig. [51 the data collapse is excellent. 

The coupling strength J affects the finite size scaling 
in two ways. First, as J increases, the interface becomes 
flatter. Accordingly, as seen in Fig. [HI for larger coupling 
strengths, the peak in fc{L) occurs for smaller targets. 
The second effect of J is to alter the depinning threshold, 
so that as J increases, /oo decreases. Heuristically, this 
can be understood by considering that as J increases, so 
does the forward force on sites lagging behind the mean 
interface position, so that the external driving force needs 
only be able to overcome the local pinning force on a few 
sites to start a large avalanche. 

For completeness, we have also made simulations us- 
ing Gaussian distributed quenched disorder, with [rji) = 
0.5,cr(?7i) = 0.2, and find very similar behavior to that 
observed with uniform-distributed disorder. 



IV. EFFECTS OF HARDENING AND 
MACHINE STIFFNESS 

Finally, we discuss two related scenarios: strain hard- 
ening and strain-controlled experiments. In materials 
subject to hardening, the external stress required to in- 
duce plastic deformation grows with the plastic strain. A 
simple model for this can be incorporated into the model 
for interface depinning, by adding a force term that is 
opposed to the interface motion, and proportional to the 
strain associated with the interface displacement, ujl? ^ 
where u is the average of Ux over x. The total force on 
site X is then 



where A; is a hardening constant. The hardening force, 
/if = —kujl?, effectively reduces the external force act- 
ing on the dislocation as the dislocation displacement in- 
creases, and for large displacements the hardening domi- 
nates and the driving force required to move the interface 
is proportional to the mean displacement. As a result, 
the interface remains critical and never depins. 

Alternatively, in the case of a strain controlled experi- 
ment, one can imagine that a subsystem undergoing plas- 
tic deformation is coupled to a spring of stiffness fc, and 
a total deformation 77^ is imposed (as illustrated in the 
Inset of Fig. [Ijb)). In the absence of a plastic subsys- 
tem, 7t would be attained using a stress of ctt — k^T, 
per Hooke's law, but when plastic deformation occurs, 
the spring extension/compression is reduced by u/Lp' so 
that ctt is given by 



(TT = A;(7T - -p). 



(6) 



Simulations of an interface subject to strain hardening 
can be rein terp reted in terms of strain controlled exper- 
iments d^, [23 , by treating / as k^T and plotting ut vs 
'jT, as illustrated in Fig. [TJb). We work in units where 
stresses and forces are interchangeable. 

We will first discuss how size effects in interface depin- 
ning are affected by strain hardening, and then discuss 
size effects in strain controlled experiments. 



A. Hardening 

We have simulated the quenched Edwards- Wilkinson 
model with an added hardening term. A typical "stress- 
strain" curve is shown in Fig. [IJa). For small driving 
force and interface displacement, the behavior is that of 
the quenched Edwards- Wilkinson model without hard- 
ening, but as the driving force and displacement are in- 
creased, the "stress-strain" curve departs from its zero- 
hardening counterpart. When hardening dominates, if a 
force of /* has been applied to drive the interface to po- 
sition u* and it has come to a stop, the force required to 
start it moving again is 



(7) 



fx = felix, Ux) - ri[x, Ux) + f - k 



L2' 



(5) 



This is the behavior seen in the "stress-strain" curve, 
which approaches a line of slope k/L^ in the limit u — > 00. 

In Fig. [71 we plot typical /c vs L curves. With the 
introduction of hardening, it is no longer possible to col- 
lapse data for the whole space of parameters 7 and k. 
Instead, the observed behavior depends on kj. For large 
k"f, fc{L) curves from simulations with the same k can 
be collapsed by subtracting k^ from (/c) , but data from 
simulations with different k cannot be collapsed onto a 
single form. Furthermore, if the k^ is small, this proce- 
dure does not yield data collapse, as exemplified in Fig. [3 
by the fc{L) — kj curve for 7 = 0.01, k = 1. 
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FIG. 7: Size dependence of depinning force for the quenched 
Edwards- Wilkinson equation subject to hardening of strength 
k. For sufficiently large kjL^ , data from simulations with dif- 
ferent 7 values can be collapsed by subtracting kj from (/c), 
to give finite size scaling of the form described by Eq. ([SI, as 
indicated by the black solid line. This collapse does not work 
for small k^L^, as seen in the data for k = l,y — 0.01 (open 
circles), which only fall on to the collapsed curve for suffi- 
ciently large L. The dashed horizontal line indicates foo, the 
L — ^ oo limit of fc{L) in the absence of hardening. Solid sym- 
bols represent A; = 10 and open symbols represent k — 1; the 
shapes represent 7 — 0.01 (circles), 7 = 0.15 (triangles) and 
7 = 2 (inverted triangles). Error bars represent the standard 
error. 

We can interpret this by considering the maximum 
hardening force applied during simulations, which is ap- 
plied when the interface is at its target: u — 7^^. For 
small ^7, the hardening force is always negligible, and 
the size effects are the same as those discussed in previ- 
ous Sections. When k"/ is large, hardening becomes im- 
portant. A first approximation to understand its role is 
to imagine that the interface is driven to the target with- 
out hardening, and then the hardening force is turned 
on. Assuming the target is sufficiently large, the force 
required to attain it (with zero hardening) is given by 
Eq. dD). Then, from Eq. ([7]), the force required to con- 
tinue to move the interface is 

f^foo+aL-^/^' + kj, (8) 

where /oo is the L ^ 00 limit of fc{L) when k = 0. 
This form agrees well with our simulation data, as seen 
in Fig. [T] Because of the extra 7-dependent term, it is 
not possible to use linear rescalings to eliminate the 7 
dependence of (/c — foo), and exact data collapse is not 
possible. 

B. Strain controlled experiments 

In strain controlled simulations, the choice of yield 
strain is more crucial than in stress controlled simula- 
tions. In Fig. ISl we plot the size dependence of (TT,y for 
several choices of yield strain jT.y, for strain controlled 
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FIG. 8: Size dependence of the measured yield stress aT,y 
in strain controlled simulations of the quenched Edwards- 
Wilkinson model, with spring constant fc = 10 (open sym- 
bols), and k — 1 (filled symbols). The yield strains used are 
7T,y = 1 (triangles), 0.5 (circles), 0.4 (inverted triangles) and 
0.04 (squares). The black lines indicate the expected finite 
size scaling of Eq. with /oo = 0.424 and 1/u = 0.75. 
When kjT.y is too small, this scaling breaks down because 
the interface never reaches criticality, as seen for the lower 
two curves {k = 10,7t,j; = 0.04, and k = l,7r,K = 0.4). 



simulations of an interface described by the quenched 
Edwards- Wilkinson equation. For sufficiently large jT,y, 
the behavior described in Eq. ^ is recovered and ax^y {L) 
approaches the same L — )• 00 limit found in stress con- 
trolled simulations, with the same finite size exponent 
l/iy — 0.75. However, for smaller jx.y, this picture breaks 
down. 

We can interpret this as follows. Recall that in stress 
controlled simulations of the quenched Edwards Wilkin- 
son model, the target -fL'^ grows faster than the interface 
width wL'^ and therefore as L — ^ 00 the correct finite size 
scaling must occur. However, in strain controlled simula- 
tions, the target "fT,y — fy/k may correspond to a driving 
force fy that is always smaller than the true depinning 
force (as given in Eq. In that case, the interface 

never reaches criticality for any L, and the standard finite 
size scaling does not hold. Note that the minimum strain 
target depends on both the depinning force of the inter- 
face and also the spring constant k, as shown in Fig. |S1 

If the target is large enough for the interface to reach 
criticality, then the size effects can be predicted by re- 
interpreting the strain controlled simulation as a stress 
controlled simulation with hardening. The finite size be- 
havior of f — kjT is given by Eq. ([1]) and ([7]), so that 

'JT,y^ f-k'-^^f 00 + aL-'/'', (9) 

where Uy is the mean interface displacement when the 
total strain is equal to the yield strain. In other words, 
for sufficiently large yield strains, stress controlled and 
strain controlled simulations give the same results, and 
furthermore, in that case, the measured yield stress in 
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a strain controlled experiment does not depend on the 
yield strain. 

V. CONCLUSION 

The question of how to measure yield stress is an im- 
portant methodological problem in the determination of 
finite size effects in microscale plasticity, which to our 
knowledge has received little attention previously. In this 
paper, we have gained insight into how to avoid spurious 
size effects in a simplified model for the onset of plastic 
yield when the yield strain is arbitrarily defined. Reas- 
suringly, in all the models we study, for sufhciently large 
systems the correct finite size scaling is recovered. How- 
ever, for smaller systems, the finite size effects are not 
described by the scaling law ([1]), and instead depend on 
various details of the model under study, a result which 
has implications especially for simulation studies of de- 
pinning and plasticity, where until recently only relatively 
small sizes were accessible computationally. 

We also observe that the idea that the crossover be- 
tween "small" and "sufficiently large" systems is deter- 
mined by the roughness of a self-affine interface is one 
that has broader relevance. For example, experimen- 



tal studies of polycrystalline copper samples under ten- 
sion |28j have revealed that profiles of the surface height 
y{x){with X along the sample axis) are self-affine, that is, 
W{L) = {\y{x + L)-y{x)\) - L^, with H « 0.75 for im- 
posed strains > 10%. This surface profile is related to the 
variation in the strain, which in a sample of length L is 
A7 = W{L)/L ~ L""-^^. Imposing the requirement that 
the total strain should be larger than its variance, a min- 
imum system size Lc is determined by setting 7 = A7. 
In the case of polycrystalline copper, this length scale 
turns out to be small enough to have no experimental ef- 
fect: using the data from Zaiser et a/.'s experiments [28} , 
Lc ~ 1.6 X 10~^ m for 7tot = 10%. For different mate- 
rials and loading conditions, however, it may be possible 
for Lc to be comparable with experimental sample sizes, 
leading to intricate size effects. 
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